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Abstract. This note deals with two topics of linear algebra. We 
give a simple and short proof of the multiplicative property of the 
determinant and provide a constructive formula for rotations. The 
derivation of the rotation matrix relies on simple matrix calcula- 
tions and thus can be presented in an elementary linear algebra 
course. We also classify all invariant subspaces of equiangular ro- 
tations in 4D. 

1. Introduction 

This article aims to promote several geometric aspects of linear alge- 
bra. The geometric motivation often leads to simple proofs in addition 
to increasing the students' interest of this subject and providing a solid 
basis. Students with a confident grasp of these ideas will encounter lit- 
tle difficulties in extending them to more abstract linear spaces. Many 
geometrical operations can be rephrased in the language of vectors and 
matrices. This includes projections, reflections and rotations, opera- 
tions which have numerous applications in engineering, physics, chem- 
istry and economic, and therefore their matrix's representation should 
be included in basic linear algebra course. 

There are basically two attitudes in teaching linear algebra. The 
abstract one which deals with formal definitions of vector spaces, lin- 
ear transformations etc. Contrary to the abstract vector spaces, the 
analytic approach deals mainly with the vector space M n and provides 
the basic concepts and proofs in these spaces. 

However, when the the analytic approach deals with the definition of 
linear transformations it uses the notion of representation of the matrix 
of the transformation in arbitrary basis and thus it actually goes back 
to the abstract setting. To clarify this issue, let us consider for example 
the calculations a rotation T of a vector x in R 3 . In this way one start 
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with choosing an appropriate basis B, and calculating the matrix of 
the given transformation in this basis [T]#. In the second step one has 
to compute the transformation matrix P from the standard basis to 
the basis B and its inverse P~ l . The final step consists of applying 
the matrix P~ 1 [T]bP to x. This cumbersome machinery is common 
for both the abstract and the analytic approaches. In addition, it take 
an essential effort to teach all the necessary details in order to use this 
non-practicable formula. 

It is amazing why one should use this complicate formula while the 
Rodrigues' rotation formula does it efficiently Although its simplic- 
ity, Rodrigues' formula does not appear in the current linear algebra 
textbooks. Its proof is elementary, but require non-trivial geometric 
insight. 

In this note we present a new proof of the matrix's representation of 
Rodrigues' formula. The essential point is that we regard the multipli- 
cation Ax, of a vector x by a matrix A, simultaneously as an algebraic 
operation and geometrical transformation (similarly to Lay [.']). This 
enable us to derive the rotation formula in the three dimensional space. 
In higher dimensional spaces we first propose a geometric definition of 
a rotation and after that we derive the formula in a similar manner to 
the three dimensional spaces. Having calculate the matrix of rotation 
according to that definition, we show it is identical to the common 
definition of rotation, that is, an orthogonal matrix with determinant 
one. 

We also consider the multiplication of two matrices AB as multipli- 
cations of the columns of B by the matrix A. Applying this point of 
view we provide a simple proof of the multiplicative property of the de- 
terminants. The standard proof of this property is often being skipped 
from the class room since it is considered as too complicated. The 
proof which we present here could easily be thought in the beginning 
of a linear algebra course. 



2. Basic facts and notations 

We recall the definitions of multiplication of a vector by a matrix 
and the multiplications of two matrices. Both definitions rely solely 
on the basic two operations of vectors in ~R n , namely, addition and 
multiplication by a scalar. 
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We denote mxn matrix A by [ai, a 2 , . . . , a n ], where {ai, a 2 , . . . , a n } 

~x{ 

are the columns of A. Let x = \ be a vector in MP, then 

_x 

(1) Ax. = XlBLl + x 2 a 2 + . . . + X n 2L n . 

This means that Ax is simply a linear combination of the columns 
of the matrix A. Conversely, any linear combination of n vectors 
{ai, a 2 , . . . , a n } can be written as a matrix multiplication. Note that 
beside of Ax being a linear combination, we can interpret it as a trans- 
formation from W 1 to M m by corresponding to each x G IR n the vectors 
Ax = y G M m . Obviously this operation has the linearity property: 

(2) A(au + /3v) = aAu + /3Av. 

Thus any linear transformation can be written as a multiplication of 
vectors by a matrix and therefore the formal definition the linear trans- 
formations from M. n to M. m seems to be superfluous. 

Most of the textbooks define a matrix multiplications by the row- 
column rule. But the original definition of Cayley is by means of a 
composition of two linear substitutions (see e.g. [4, 6]). This means 
that if A is a m x n matrix and B n x k, then the matrix C = AB is 
defined through the identity Cx = A(Bx), where x G M. k . From (1) 
and (2) it immediately follows that 

(3) AB=[Ah 1 ,Ab 2 ,...,Ah k ], 

where {b l5 b 2 , . . . , b^} are the columns of the matrix B. 

Another useful way to multiply matrices is by the column row rule, 
that is, 

(4) AB = a x b^ + . . . + a n b^, 

where , . . . ,b^ are rows vectors of the transpose B T . This type of 
multiplication will be used in Section 4. 



3. Multiplicative property of determinant 

Let Abe nx n matrix with coefficients [«ij]^ =1 - The determinant of 
A is defined to be the scalar 

(5) det(A) = 2^ o-(p)a lpi a 2p2 . . . a npn , 

p 
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where the sum is taken over the n\ permutations p = (pi,p 2 , • • -Pn) of 
(1,2, ... ,n) and 



a(p) 



-1 if pis even permutation 
-1 if pis odd permutation 

Theorem 1. Suppose A and [u, v, . . . , z] are two n x n matrices, then 
(6) det ([An, Av, . . . , Az]) = det(A) det([u, v, . . . , z]). 

Proof Denote A by [a l5 a 2 , . . . , a n ] and let 
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We are now using formula (1), the linearity of the determinant and the 
fact that a matrix having two equal columns its determinant is zero. 
All these result with 



det ([Ail, Av, . . . , Az]) = det 



i=l 



1=1 



i=l 



— ^^(u Pl v P2 . . . z Pn ) det([a Pl , a P2 , . . . , a Pn ]). 

p 

Since the determinant changes sign when two columns are interchanged, 
we have 



det ([a 



Pi 1 ™~P2 1 



•••, a P„]) = ^(p) det (A). 



Therefore, 



U Pl V P2 



-Vn 



det([Au, Av,..., Az]) = det(A) a(p) 

p 

= det(A) det([u, v, . . . ,z]). 

□ 

The above Theorem has an important geometric interpretation and 
we discuss it here in M. 2 . It is well known that if u, v are two collinear 
vectors, then det([u, v]) is the area of parallelogram spanned by {u, v}. 
Now if det (A) ^ 0, then {Au, Av} span also a parallelogram. There- 
fore the number det (A) is the proportion between the areas of the 
parallelograms spanned by {u, v} and -{Au, Av}. 

Note that Theorem 1 gives the multiplicative property of the deter- 
minant. Indeed, let B — [u, v, . . . , z], then AB = [Au, Av, . . . , Az] and 
(6) becomes 

det (AB) = det (A) det(B). 
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4. The derivation of a formula for 
a rotation matrix in w 1 

First of all let us discuss the definition of a rotation in R n . The com- 
mon definition of a rotation is by means of an orthogonal matrix with 
determinant one. Here we provide another definition based on geo- 
metric considerations. We are aware that such definition was probably 
given in the past, but we could note traced it. Its advantage is being 
practicable and the computation of the rotation's matrix does not uti- 
lize eigenvalues and eigenvectors. We shall then verify the equivalence 
of the two definitions. 

It is rather simple to define a rotation in IR 2 . A linear transformation 
R is a rotation if the angle between the vectors i?x and x is a constant 
for all x G M 2 . From this definition follows that 

R 

where the rotation is by an angle a counterclockwise. Therefore 
(7) 



"1 


0~ 




cos a 


— sin a 





1 




sin a 


cos a 



R 



cos a — sin a 
sin a cos a 



cos al + sin a J 



where I is the identity matrix and J 







Formula (7) implies 



that confirms our geometric 



that ||-Rx|| = ||x|| for each x G IR and 
intuition. Matrices which preserve norm are called orthogonal matrices. 
Their determinant is ±1 and rotations are orthogonal matrices with 
determinant one. 

The definition of a rotation in IR 3 is slightly more involved. A linear 
transformation R is called rotation if there exists two dimensional sub- 
space II of IR 3 (the plane of the rotation) such that the angle between 
vectors i?x and x is a constant for all x G II, and Ry = y for each y 
orthogonal to IT (the axis of the rotation). Euler's theorem about the 
rigid motion of the sphere with fixed center justifies this definition. 

In order to calculate the matrix R we pick two orthonormal vectors 
a, b G IT and a unit vector c which is orthogonal to II and such that 
the triple {a, b, c} is right-handed. Applying the rotation R to these 
vectors results with 



R& = cosaa + sinab, Rh 



sin aa + cos ab, Rc = c. 



We write the above equalities in a matrix form RP = Q, where P = 
[a, b, c] and Q = [R&, Rb, Rc}. Since P is an orthogonal matrix, R = 
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QP T and calculating QP T by means of (4), we get that 



R = [Ra, Rh, Rc] 



a T ' 



{Ra)a T + (Rb)b T + (Rc)c T 



= cos a (aa T + bb T ) + sin a (ba T — ab T ) + cc T . 

The skew symmetric matrix (ba T — ab T ) is the matrix representation 
of the cross product c x x. To see this note that (ba T — ab T )c = 0, 
(ba T — ab T )a = b and (ba r — ab T )b = —a. Hence 

i?x = cos ax + (1 — cosa)cc T x + sina(c x x), 

which is the known Rodrigues' formula. 

Formula (8) resembles to a large extent the two dimensional formula 
(7). The matrix (aa T + bb T ) is the projection on the plane II, cc T is 
the projection on the line orthogonal to II and 

(9) (ba T - ab T ) 2 = -(aa T + bb T ). 

Since the rotation is actually in the plane IT, we see that (aa T + bb T ) 
corresponds / and (ba r — ab T ) corresponds J in formula (7). In both 
formulas, cos a is the coefficient of a symmetric matrix and sin a is the 
coefficient of an anti-symmetric matrix. 

It is easy to check that the matrix R(a) := R in (8) is an orthogonal 
matrix with determinate one. Indeed, relation (9) and the orthogonal- 
ity of {a, b, c} imply that 

R(a)R T (a) = (cos 2 a + sin 2 a) (aa T + bb T ) + cc T = I 

and hence det(R(a)) = ±1. Letting lim Q ^ R{oc) = I and using the 
continuity of the determinants, we see that det(R(a)) = 1. 

We turn now to rotations in M 4 . It turns out that rotations in R 4 
can be defined in a similar way to rotations in M 2 and M. 3 . We say that 
a linear transformation R is a rotation if there exists two dimensional 
subspace II of IR 4 such that the angle between vectors i?x and x is a 
constant for all x G II, and the angle between vectors Ry and y is a 
constant for all y G H ± , the orthogonal complement of II. 

The calculation are done in a similar manner as we did in R 3 . Pick 
a, b G II and c, d G H ± such that the set {a, b, c, d} is an orthonormal 
basis. Let R = R(a, (3) be the rotation matrix with rotation's angles a 
in the plane II and f3 in the orthogonal complement il- 1 . Then 

Ra = cos aa + sin ab, Kb = — sin aa + cos ab, 
Rc = cos /3c + sin (3d, Rd = — sin j3c + cos j3d. 
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Set P = [a, b, c, d] and Q = [Ra, Rb, Rc, Rd], since the matrix P is 
orthogonal, R = QP T and hence 

R = (i?a)a T + (Rb)b T + (Rc)c T + (Rd)d T 
(10) = cos a (aa T + bb T ) + sin a (ba T — ab T ) 

+ cos /3 (cc T + dd T ) + sin /3 (dc T - cd r ) . 

We can now easily distinguish between two types of 4D-rotations. If 
(3 = 0, then the rotation is simple, that is, Ry = y for all y G il -1 . 
Otherwise, both planes II and II 1 - rotate simultaneously and this type 
is called a double rotation. 

If one doubts whether the matrix R in (10) is an orthogonal ma- 
trix with determinant one, then the following simple calculation will 
convince him. Since 

R(a,P)R T (a,/3) 

= (cos 2 a + sin 2 a) (aa T + bb T ) + (cos 2 /3 + sin 2 fj) (cc T + dd T ) 
=/, 

R(a, j3) is an orthogonal matrix and by letting a and /? go to zero, we 
get that its determinant is one. 

We are now in position to extent the geometric definition of rotations 
to M n for arbitrary positive integer n. For n = 2p we say that a linear 
transformation R is a rotation if there exist p mutual orthogonal planes 
Uk such that the angle between the vectors i?x and x is a constant for 
all x G Uk, k = 1, ...,p. For n = 2p + 1 we require that there are p 
mutual orthogonal planes LT^ and in addition a line L is orthogonal to 
Uk, k = 1, . . . ,p such that R behaves the same as in the even on the 
planes 11^ and Ry = y for all y G L. The extension of formulas (8) and 

(10) to arbitrary dimension is obvious. In M 2p there is an orthonormal 
basis {(ai, bi), (a n , b n )} such that 

p 

(11) R = 22 cosa k (a fc a^ + b k bl) + sma k (b k &l - a fc b^) . 

k=l 

In odd dimension 2p + 1 there is an orthonormal basis {(a l5 bi), 
(a„, b n ), c} such that 

p 

(12) R = ^ cos a k ( a fe a fc + b fcbfc ) + sin a k (b fc a^ - a fc bjQ + cc T . 

k=l 

Similarly to the rotation formulas (8) and (10) one can check that 
(11) and (12) are orthogonal matrices with determinate one. 
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Formulas (11) and (12) were derived by [5] but in a different way. 
The advantage of the derivation given here is being constructive an 
addition to being appropriate for an elementary linear algebra course. 
Formulas (11) and (12) can be written in a vectors' form 

v 

Rx = ^2 cos a kYk + sin a k z k 
fc=i 

and 

p 

#x = ^ cos a k y k + sin a k z k + (c T x)c, 

fe=i 

where y& is the projection of vector x on the plane 11^ and z k is the 
rotation of y^ by an angle | in the plane 11^. 

4.1. Invariant subspaces of equiangular subspaces of rotations 

in R 4 . A rotation R in IR 4 is called equiangular rotations or isoclinic 
rotations if the planes II and its orthogonal complement Il ± rotate with 
the same angle (see e.g. [2, 3]). When a ^ (3, then the planes II and U 1 - 
are the only invariant subspaces under the rotation R. However, when 
a — /3, then there are infinitely many two dimensional invariant planes 
(see e.g. [5]). We shall see here that this interesting phenomenon is 
a simple consequence of the formula (10) and we shall also classify all 
the invariant planes. 

To see this we note that when a — j3, then (10) becomes 

R = cos al + sin a J, 

where / is the identity on IR 4 and J = (h& T — ah T + dc T — cd T ) . 
Hence a subspace U of IR 4 is an invariant subspace of R if and only if 
it is invariant subspace of the matrix J. 

Now J is a skew-symmetric matrix satisfying J 2 = —I. Therefore it 
has no real eigenvalues and this implies that any non-trivial invariant 
subspace must has dimension two. Since J 2 = — /, span{u, Ju} is an 
invariant subspace for any u e IR 4 . On the other hand, if U is a non- 
trivial invariant subspace of R, then Ju £ U for any u £ U. Hence U 
must be spanned by these vectors. Thus we have obtained a complete 
classification of the invariant subspaces of equiangular rotations which 
is independent of the rotation angle a. 

It follows from formulas (11) and (12) that if all the angles a k are 
equal, then there are infinitely many invariant subspaces and each one 
of them is spanned by a vector u and Yl k =i (bfc a fc — a fcbfc ) u. 
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5. Concluding Remarks 

Many vector-space textbooks use the entry-by-entry definition = 
a ikbkj for the matrices multiplications. The operation of multiplica- 
tion of a vector x by a matrix A in accordance (1) bears in itself both 
geometric and algebraic properties. Therefore a decent understanding 
of it should be prior to the formal definition of matrices multiplication. 
After that the matrices' multiplication in Cayles's spirit follows nat- 
urally. The Cayley's definition (3) and the column-row rule (4) have 
many advantages. In many cases they makes the computations easier 
in addition to increases the comprehension. This note emphasizes two 
aspects of that attitude. 
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